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TITLE; 

TYPE; 

LEK3TH: 

TEMPORARY STORAGE; 



Romberg integration 

closed^ relocatable^ mnemonic 

2$ words 

6 words at fixed memory locations 0-5 
n words beginning at location S 



FAST REGISTERS CHARGED; none 
SUBROUTINES TEED; none 
DURATIOH: 



PARAMETERS z 



depends on the duration of the auxiliary routine 
for evaluating the function f(x). Since the 
auxiliary is entered 2 ~ + 1 times, duration is 
approximately doubled when n is increased by 1. 



link in ML5 

k parameters which must be written in the word 

following the word containing the JSB instructions; 



2 < n < 13 



a 




n 



address of left limit a of integration interval 

address of auxiliary routine 

address of first location of temporary storage 

v\ 1 

integration interval divided into 2 
subintervals 



SESCRXPTXOls 



The right limit b of the integration interval has to 
be placed in the location immediately follo^ring the 
one which contains the left limit a. 

this subroutine computes I = / f(x) ax 
in the following way; 



-1- 



Compute' a sequence R f R . , ... R of approximations to I, where R is 
found by dividing the integration interval (a,b) into 2 " sub intervals of 
equal length and evaluating the integral by the trapezoidal rule, 
The sequence R.„ converges linearly to I with convergence factor 1/4, i.e, 

H i,l " X - 1/k 

i 1 " K i-1 1 
■Ehe new sequence R^ = _*_^- 

also converges linearly to I, with convergence factor l/l6„ 



Beginning with R . we generate different sequences R which all converge 
linearly to I according to the formula 

and arrange them as columns in a triangles 
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The user must provides 

The left limit a and the right limit b of the integration interval^, 
la tm> consecutive locations^ and specify the address of a 
as parameter j, 

An auxiliary routine for evaluating the function f, ishieh is also 
linked in ML5, 

A block of n consecutive words in msmory "beginning at location S. 

The subroutine stores R^,/("b»a) B^/Cb-a) ... B^/Cb-a) 

in memory locations S S+l S+n-l 

and leaves the best approximation H to the integral in the 

nn 

accumulator. 



The choice of the parameter a {'which implies that the integration 
interval is divided into 2 "' subintervals) should depend on the 
behavior of the function f(x), taking into consideration that 
the duration is approximately doubled vhen n is increased by 1. 



Range; 2 < n < 13° 



BATE: December l s 196*2 

BY: J. RUvergelt 
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read parameters Into F8 
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